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ABSTRACT 

The evidence for a shortage of exosolar planets with semimajor axes — 1.1 < 
\og(a/AU) < —0.2 is investigated. It is shown that this valley results from a gap 
in the radial distribution of planets, orbiting stars with masses > 1.2M (the 
high-mass sample, HMS). No underabundance is found for planets orbiting stars 
with smaller masses. The observational data also indicate that within the HMS 
population it is preferentially the more massive planets with Msin(z) > 0.8Mj 
that are missing. Monte-Carlo simulations of planet formation and migration are 
presented that reproduce the observed shortage of planets in the observed radius 
regime. A dependence on the disk depletion timescale Td ep is found. The gap is 
more pronounced for Td cp = 10 6 — 10 7 yrs than for Td ep = 3 x 10 6 — 3 x 10 7 yrs. 
This might explain the observed trend with stellar mass if disks around stars 
with masses M* > 1.2M have shorter depletion timescales than those around 
less massive stars. Possible reasons for such a dependence are a decrease of disk 
size and an increase of stellar EUV flux with stellar mass. 

Subject headings: planetary systems: formation - extrasolar planets 

1. Introduction 

Within the past decade, nearly 200 extrasolar planets have been detected. Their orbital 
properties and masses and the relationships between planetary properties and those of the 
central stars are providing interesting insight into the complexity of planetary system for- 
mation and offer important constraints for theoretical models. Recently, extrasolar planets 
have been found around low-mass stars (M dwarfs) by high-precision radial-velocity (Butler 
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et al. 2004, Bonfils et al. 2005) and microlensing surveys (Beaulieu et al. 2006, Gould et 
al. 2006) with indications that Jupiter-mass planets are rare around M dwarfs (Laughlin et 
al. 2004, Ida & Lin 2005). 

Theoretical models have been developed to explain the statistical properties of exoplan- 
ets. Ida & Lin (2005) presented Monte-Carlo simulations that predict mass and semimajor 
axis distributions of extrasolar planets around stars with various masses. Their prescrip- 
tion of planet formation (Ida & Lin 2004a) is based on the core-accretion scenario in which 
it is assumed that Jupiter- mass gas-giant planets form through 1) grain condensation into 
kilometer-sized planetesimals, 2) runaway planetesimal coagulation (Wetherill & Stewart 
1989; Kokubo & Ida 1996), 3) oligarchic growth of protoplanetary embryos (Kokubo & Ida 
1998, 2000), and finally 4) gas accretion onto solid cores (Mizuno 1980; Bodenheimer & Pol- 
lack 1986; Pollack et al. 1996; Ikoma et al. 2000). The models show that the characteristic 
mass distribution of planets around M dwarfs should be quite different from that around 
FGK dwarfs, which agrees well with the results of microlensing surveys (Beaulieu et al. 2006; 
Gould et al. 2006). The core accretion model also predicts a strong dependence of the like- 
lihood for finding planets on the stellar metallicity (Ida & Lin 2004b; Kornet et al. 2005) 
which is in agreement with observations (Santos et al. 2000; Fischer & Valenti 2005). 

Like the mass distribution, the semimajor axis distribution of exoplanets needs to be 
explained by any consistent theoretical model of planet formation. One of the most puzzling 
questions that arose already when the first Jovian planets around main sequence stars were 
discovered (Mayor & Queloz 1995) is their wide range of orbital radii, with distances ranging 
from 2 x 10 -2 AU to the observational limit of ~ 5 AU. According to the core-accretion 
scenario, preferred formation locations of Jovian planets are the outer regions of protoplan- 
etary disks where the temperatures are low enough for ice to survive (e.g. Ida & Lin 2004a). 
However, hot Jupiters have been observed very close to the central star (< 0.1AU). It is 
unlikely that these objects formed in these inner regions, although Jovian planets can form 
inside the ice boundary in relatively massive disks (Ida & Lin 2004a). Orbital migration due 
to gravitational interaction between the planet and the gaseous disk (Goldreich & Tremaine 
1980; Lin & Papaloizou 1986a,b; Lin et al. 1996) is considered as the most promising mech- 
anism to generate hot Jupiters. It starts with relatively fast type I migration (e.g. Ward 
1986; Tanaka et al. 2002). As soon as the planet has become massive enough to open a gap 
in the disk it enters the phase of type II migration where the infall rate is connected with the 
viscous evolution of the disk. Numerical simulations have explored the complex migration 
process in great details (see e.g. de Val-Borro et al. 2006 for a recent comparison of different 
numerical simulations) . 

The Monte Carlo models of Ida & Lin (2004b) predict that the (differential) distribution 
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for extrasolar planets with radial velocities larger than 10 m/s is an increasing linear function 
in log a from a ~ 0.1 AU up to ~ 1 AU and decreases with a for a > 1AU, independent of disk 
metallicity, if the disk depletion timescale is comparable to the disk diffusion timescale at the 
radius where Jovian planets are formed (The diffusion timescale determines the migration 
speed). Very similar results were found by Armitage et al. (2002) and Trilling et al. (2002). 
Some planets might survive in the outer regions. Near the inner disk edge (< 0.1AU), those 
planets that migrate inwards from the outer regions may pile up (Lin et al. 1996). As a result, 
a gap (valley) would be created at intermediate semimajor axes (equivalently, intermediate 
periods) in the distribution. Udry et al. (2003) indeed point out a shortage of planets in 
the 10-100 day period range in the observed data (see also Zucker & Mazeh 2002) which 
they interpret as a signature for a transition region between two categories of planets that 
suffered different migration scenarios. 

In this paper we analyse again the valley in the period and separation distribution 
of extrasolar planets taking into account the newest set of observational data. We point 
out that in the observational data, the valley is much more pronounced around F dwarfs 
than GK dwarfs which may reflect a dependence on disk size and/or on the EUV flux from 
the central star. Section 2 presents the observational evidence for a gap around F dwarfs. 
Section 3 discusses a series of Monte-Carlo simulations that lead to a dependence of the 
planet semimajor axis distribution on stellar mass. A discussion of the results follows in 
section 4. 



2. A gap in semiaxis distribution for planets around higher-mass stars 

We use the sample of extrasolar planets around normal stars, compiled by the Extrasolar 
Planets Encyclopaedia (http://exoplanet.eu/) and include all planets (total number: 175) 
with known semimajor axes a and stellar masses M* that have been detected by radial 
velocity measurements. The sample is separated into 8 bins, equally spaced in log a between 
-1.7 < log (a/ AU) < 0.7. 

We now split the sample into planets (38) orbiting stars with masses M* > 1.2M Q 
(which we denote the high-mass sample, HMS) and the rest (137) which orbit less mas- 
sive stars (the so called low-mass sample, LMS). The right panel of figure 1 shows the 
distribution of semimajor axes of the LMS. The planets are homogeneously distributed for 
— 1.7 < log(a/AU) < —0.2 with a larger frequency in the three bins with large semimajor 
axes, corresponding to log (a/AU) > —0.2. No clear valley is evident for log (a/ AU) < —0.2. 
The distribution looks however very different for the HMS (left panel of figure 1) where a 
clear lack of planets is found in the radius regime of —1.1 < log (a/AU) < —0.2. 
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We tested the assumption that the planets in the 5 inner logarithmic bins with a < 
0.6 AU are consistent with a homogeneous distribution. The chi-square test of the LMS 
planets leads to a likelihood of 78% that this sample is homogeneously distributed in log a. 
In contrast, the chi-square statistic of the HMS planets gives a likelihood of only 2% that 
the gap in the semimajor axis distribution is a result of statistical fluctuations. Although 
the statistical significance is still marginal due to the low number statistics, there is a good 
chance that the lack of planets for stars with masses above 1.2 M in the semimajor axis 
regime 0.1 AU < a < 0.65 AU is real, while the distribution is much more uniform in log a 
for less massive stars. 

The lower two panels of figure 1 show an even more detailed analyses of the planets in 
the HMS. In the left panel all planets with masses Mshu > 0.8Mj are plotted. A large gap 
is visible. In contrast, the distribution of planets with Msini < 0.8Mj is homogeneously 
distributed with no gap (lower right panel of figure 1). Note also in the lower right panel a 
lack of planets with large semimajor axes log (a/AU) > —0.2, which was also discussed by 
Udry et al. (2003). 



One of the possibilities to account for the different planetary semimajor axis distribu- 
tions of the HMS and LMS is different sizes of their protoplanetary disks. Armitage et 
al. (2002), Trilling et al. (2002) and Ida & Lin (2004b) find that the semimajor axis dis- 
tributions are controlled by the ratio of the disk depletion timescale (7d cp ) to the viscous 
diffusion timescale T v (ri), evaluated at the formation site of the Jovian planets {r\ ~ 1- 
10AU). Because T v {r\) determines the type II migration speed of Jovian planets while Td ep 
determines the duration of type II migration, Jovian planets migrate more and the semima- 
jor axis distribution is smoothed out more for larger Td cp /r I/ (ri). For smaller rd ep /r I/ (ri), on 
the other hand, Jovian planets tend to retain their original locations and the semimajor axis 
distribution shows more clear gaps. We show below that r dcp /r I/ (r 1 ) decreases with stellar 
mass, leading to a more pronounced gap for higher-mass stars. 

The formation of Jovian planets is favored in the regions just beyond the ice boundary 
(e.g., Ida & Lin 2004a). It is therefore reasonable to identify r 1 with the radius r icc of the 
ice boundary (Hayashi 1981) 



where we assume that the stellar luminosity is proportional to M*, although the depen- 
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dence may be slightly weaker for pre-main sequence stars. We adopt a disk temperature T 
that is regulated by stellar irradiation (Hayashi 1981), 

T = 280 1/2 (^ K. (2) 

Then the viscosity v = ac 2 s /Vt oc arM^ 2 , where c s oc \[T is the sound velocity, Vt oc r~ 3 / 2 
is the Keplerian frequency and a is the parameter of the alpha-viscosity model (Shakura & 
Sunyaev 1973). With these assumptions for r x and v, 

The estimate the disk depletion timescale we adopt the standard similarity solution of vis- 
cously evolving disks (e.g., Lynden-Bell & Pringle 1974, Hartmann et al. 1998). The disk 
depletion time is then 

r dcp ~ r,(r ) ~4t~3x10 6 (-2L)" 1 f ^V^yrs. (4) 

p v o; 3i/(r ) V10" 3 / V100AU/ V^©/ 

where the disk size r corresponds to the radius of the maximum viscous coupling where the 
radial dependence of the disk gas surface density (E g ) changes from r' 1 to exp(— r/r ). For 
t < T„(ro), r is constant with time. 

Recently, the dependence of protoplanetary disks on their central stars is being obser- 
vationally explored. One of the most remarkable dependences that has been found is that 
the mass accretion rate onto the stars is proportional to the square of the stellar mass M* 
(Muzerolle et al. 2005; Mohanty et al. 2005; Natta et al. 2006). Alexander & Armitage 
(2006) point out that a decrease in disk size r with increasing (r oc M* ) can ex- 
plain such a strong dependence. In this case, Td ep oc M~ x (Eq. 4). From this and Eq. (3), 
Tde P /Tz/(^i) oc M* 5 ^ 2 . Therefore, Tdcp/V^i) may differ by a factor 3 between disks around 
stars with ~ O.9M and those with ~ 1.4M . 

Motivated by the above arguments, we have carried out Monte-Carlo simulations with 
fixed a but with various Td cp to determine the semimajor axis distribution of planets with 
radial velocities larger than 5 m/s. The details of the calculations of planetesimal accretion, 
gas accretion onto cores, gap opening and migration are described in Ida & Lin (2004a, b, 
2005). Here, type II migration with a = 10~ 3 is included, but type I migration is neglected. 
The effect of type I migration will be presented in a separate paper (Ida & Lin 2006). 
We assume that the surface density of planetesimals is Ed = /disk£d,MMSN ^ r~ 3 ^ 2 , where 
Sd,MMSN is that of the minimum-mass solar nebula model (Hayashi 1981) and the scaling 
factor /disk takes values of 0.1-10 (see discussion in Ida & Lin 2004a). On the other hand, 
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the gas surface density is S g = / disk (r/lAf/) 1 / 2 S gi MMSN(l + t/r dep ) 3/2 oc r 1 (1 +t/r dcp ) 3 / 2 , 
corresponding to the similarity solution at r < r . 

As long as FGK dwarfs are considered, the formation mechanism of planets is almost 
independent of M* and stellar (disk) metallicity (Ida & Lin 2004b, 2005). We therefore show 
results only for M* = 1M and [Fe/H] = 0.2 to make the discussion simple. With M* = 1M 
and a = 10~ 3 , the viscous timescale in all simulations is T v {ri) ~ 10 5 yrs (Eq. 3). Since the 
key parameter to determine the radial variation in the planet semimajor axis distribution is 
7dep/7v( r i), we did calculations with Td ep = 10 6 -10 7 yrs, 3 x 10 6 -3 x 10 7 yrs, and 10 7 -10 8 yrs, 
assuming that distributions of Td ep are log uniform in these ranges. Note that the results 
with shorter Td ep correspond to higher M*. For a = 10~ 2 , 10 times shorter r dcp produce 
similar results. 

The upper panels in Fig. 2 show the theoretically predicted mass and semimajor axis 
distributions. The left panel assumes r dcp = 10 6 -10 7 yrs, while the right one corresponds to 
Tdep = 3 x 10 6 -3 x 10 7 yrs. We identify the former and the latter results with planets forming 
around HMS and LMS, respectively. The middle panels are histograms of the semimajor 
axis distributions. Because we artificially terminate type II migration at 0.04AU (Ida & Lin 
2004a), we obtain too much pile- up of hot Jupiters. In reality, many of these hot Jupiters 
may actually either be consumed (e.g., Sandquist et al. 1998) or tidally disrupted (e.g., 
Trilling et al. 1998; Gu et al. 2003) by their host stars. The gap at 0.1-1AU is pronounced 
in the result of r dcp = 10 6 -10 7 yrs, while that of Td cp = 3 x 10 6 -3 x 10 7 yrs shows a flatter 
distribution. Because the decline of type II migration is earlier in the shorter Td ep case, 
initial formation locations are more frequently retained in the results with Td ep = 10 6 -10 7 
yrs than in the case of r dcp = 3 x 10 6 -3 x 10 7 yrs. In the lower panels, we split the histogram 
corresponding to r dcp = 10 6 -10 7 yrs (middle left panel) into two parts: M > 0.8Mj (lower 
left panel) and M < 0.8Mj (lower right panel). Because type II migration is slower for 
more massive planets, the gap is more pronounced for M > 0.8Mj. Thus, if the results with 
Tdep — 10 6 -10 7 years and 3 x 10 6 -3 x 10 7 yrs represent HMS and LMS, respectively, the 
theoretical predictions show a trend with stellar and planetary mass that is consistent with 
the analysis of the observational data in section 2. 

Although statistical fluctuations may still be large in the observations presented in 
Fig. 1, the rise of the mass distribution at ~ 1 AU appears to be steeper than the theoretical 
prediction. This might be due to photo-evaporation of disk gas due to EUV radiation from 
the central star which is very effective at a few AU and could inhibit type II migration for 
planets that formed beyond a few AU (Matsuyama et al. 2003). Note however that the 
process of disk gap formation by photoevaporation depends critically on the stellar EUV 
flux in the pre-main sequence phase which is poorly known. The X-ray flux from pre-main 
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sequence stars, on the other hand, is observed and the flux may be linearly proportional 
to M* (Preibisch et al. 2005). Hence, the EUV flux may also increase with M*. Although 
the dependence with stellar mass may not be strong enough, photo-evaporation effects could 
make the radial gradient of the planet mass distribution at ~ 1AU steeper for the HMS 
sample than for the LMS sample. More theoretical work would be required to understand 
this process in greater details. 

4. Discussion 

We should start the discussion with a word of caution. Despite the fact that almost 200 
exoplanets are known, the statistics for more advanced correlations like the one discussed 
here is still poor. In addition, for small numbers spurious correlations between different 
variables can always be found. They would appear to have a high statistical significance 
which is however caused by the biased selection of the data. At the moment we cannot rule 
out that the planet shortage in the separation range of 0.08 AU < a < 0.6 AU, first discussed 
by Udry et al. (2003), is such a case. 

However, if this valley in period and separation distribution exists, we argue that it is due 
to a lack of planets with masses M sin(i) > 0.8 Mj orbiting around stars with masses M* > 
1.2 M which provides interesting new insight on the origin of massive planets orbiting close 
to their parent stars and the dependence on stellar properties. We point out that smaller size 
protoplanetary disks around more massive stars can account for a deeper valley around these 
stars. If MRI turbulence is responsible for the disk viscosity, a values for the viscosity would 
not significantly depend on stellar mass. Adopting a constant a, smaller disk sizes around 
more massive stars would lead to shorter disk depletion timescales, which prevents Jovian 
planets formed in the outer regions from migrating considerably. Strong EUV radiation 
from massive (pre-main) stars may also inhibit the migration across the regions at a few 
AU. Another possibility is a larger radius of the ice boundary for massive stars that may 
also inhibit smoothing-out the valley. The valley would then be most pronounced for planets 
around massive stars. 

To confirm the existence of the HMS period valley and test the scenario presented here, 
more detections of extrasolar planets around more massive stars (BA dwarfs) are needed. If 
our model is correct, the valley in the period/separation distribution of Jovian planets would 
become even more pronounced for even more massive stars. Because several planets around 
KG giants, which were BA dwarfs in their main sequence phase, have been discovered, plenty 
of planets should also exist around BA dwarfs. Spectroscopic observations are not easy for 
BA dwarfs because of a smaller number of absorption lines and rapid rotation. However, in 
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spite of the difficulty, relatively massive planets should be detectable. Transit surveys may 
not have serious problems to search for planets around BA dwarfs, however, detection is 
limited to short-period planets. Detailed observations of differences in disk sizes between T 
Tauri stars and Herbig Ae/Be stars are also needed to test our model. 
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Fig. 1. — The semimajor axis distribution of extrasolar planets orbiting stars with masses 
M* > 1.2M (upper left panel) and stars with masses M* < 1.2M (upper right panel). The 
lower two panels show the distribution of extrasolar planets with masses M > 0.8Mj (lower 
left panel) and M < 0.8Mj (lower right panel) for the HMS sample (upper left panel). 
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Fig. 2. — The mass distributions of extrasolar planets predicted by the theoretical model. 
Upper left and right panels show the mass and semimajor axis distribution with T^ cp = 10 6 - 
10 7 yrs and 3x 10 6 -3x 10 7 yrs, respectively, which may correspond to HMS and LMS samples. 
For other parameters, see text. Middle panels are histograms of mass distributions of planets 
with radial velocity larger than 5m/s (corresponding to currently observable planets). In the 
right panel, the result with 10 7 -10 8 yrs is also plotted with open circles. Lower panels are 
histograms for planets with M > 0.8Mj (lower left panel) and M < 0.8Mj (lower right 
panel) for the r dcp = 10 6 -10 7 yrs result (middle left panel). 



